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PROBLEM 


Investigate and report on factors in physical oceanography which pertain to 


underwater sound, including physical properties such as sound velocity, thermal 
structure, and water motion. Specifically, determine the type of shallow-water 
flowmeter best suited for measuring low level, steady-state currents in the 
presence of orbital velocity motion due to surface waves. 


RESULT 


No known commercial meters are available for accurate measurements of 


this type. The electromagnetic method is believed to be inherently the best 
approach to the problem. 


RECOMMENDATIONS 


1. Investigate new configurations using ac fields that will 


(a) provide minimum hydrodynamic interference, 


(b) measure flow over large distances. 


2. Experimentally determine the directional response of the ‘‘tube”’ configuration. 


3. Consider using the earth’s field type flowmeter, but only if a feature for zeroing 


is employed. 
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INTRODUCTION 


Flowmeters are not normally designed for the purpose of measuring shallow- 
water, steady-state currents in the presence of orbital velocities. A need therefore 
exists for information on the value of existing meters used for this purpose. 

This report presents a feasibility study of the relative merits of various 
instruments and techniques. The dynamic response, filtering, linearity, and direc- 
tionality of various devices were investigated to identify the most accurate flow- 
meter to employ at depths where the magnitude of orbital velocities caused by 
surface waves can be several times larger than steady-state flows. 


PROCEDURES 


Selection 


Typical magnitudes of orbital velocities are from 5 to 50 centimeters per 
second (cm/sec) with periods from 10 to 15 seconds. Typical steady-state signals 
are actually very-low-frequency velocity fluctuations. Typical magnitudes are 
from 2 to 10 cm/sec with periods of interest at 5 to 10 minutes and at semidiurnal 
and diurnal regions. Since the frequency of these oscillations is so low compared 
with the predominant orbital or swell frequencies, they may be considered de or 
steady-state signals for most analytical work. 

The primary concern will be the selection of a transducer with inherently 
good features for this particular type of measurement. Of secondary concern will 
be items such as the amount of electronic equipment required or the necessary 
power and maintenance. The accessibility of flowmeters installed in shallow 
water (compared, for example, with deep-water buoy installations) allows more 
flexibility in the specifications. If sizable amounts of electronic equipment need- 
ing ac power are required, this will not be considered an inherent bad feature of 
the system. However, if the electronic package has to be so close to the sensor 
that it interferes hydrodynamically with the flow, then it would become an inherent 
fault. 


Operation 


Two basic procedures, continuous and discrete, can be used operationally 
for the measurement of the velocities of currents. The continuous method can use 
a transducer: 


1. With frequency response from zero to above the swell frequency; all 
data are recorded and the swell data removed later. 


2. As in 1., but with swell frequencies filtered out between transducer 
and recorder. 


3. That filters out the swell frequency and records all other data (rejec- 
tion filter). 


4. That passes only signals with frequencies lower than that at the swell 
(low pass filter). 


The best approaches would be 3. and 4. because the filtering is done at 
the earliest possible stage and a roughly flat spectrum is recorded. The choice of 
upper frequency cutoff is based on highest frequency of interest and highest fre- 
quency recorded.! A typical ratio is 2:1. There would be little difference in data 
recorded in 3. and 4. because of the low energy content of the spectrum above the 
swell frequency. Approach 1. is undesirable because of the small signal-to-noise 
ratio and large amount of data to record and reduce. 

The discrete method can be imposed (either after recording or in real time) 
or inherent. Inherent discreteness is found by timing particles traveling over 
known distances. Care must be exercised in choosing the sampling interval so 
that aliasing does not occur. 


REQUIREMENTS 


The characteristics and conditions for a successful flowmeter are as 
follows: 


1. Outputs with orthogonal components, or a capacity to measure the 
velocity directly with outputs for magnitude and angle; resultants in the horizontal 
plane are sufficient, but three-dimensional types are more desirable. 


2. Measurement errors of the steady-state value less than 5 percent. 


3. A low-temperature coefficient needed because many of the signals of 
interest are associated with large temperature changes — 5 to 10°C are typical. 


4. Meter mounts affixed to the sea bottom or stationary structures to pre- 
vent motion relative to the earth; errors resulting from mounting to buoys are not 
considered. 2 


5. Large dynamic velocity range. 


6. Bidirectional response, the angular response preferably a cosine 
function. 


7. Linear output. 


a. De levels not connected with swell amplitudes. Because non- 
linearities are a function of amplitude, the swell magnitude 
can generate different, false dc signals in some cases. This 
is particularly important when looking at diurnal signals. 


b. Intermodulation distortion from nonlinearities, where more than 
one strong frequency is present in the natural spectrum, should 
be low. 


8. Long-term stability. The ability to check and adjust (or record) the 
zero and sensitivity is a necessary feature. 


9. Electrical output. 


10. Size that will permit testing in tow tanks. 


FLOW METERS 


Acoustic 


Acoustic flowmeters are of two basic kinds: (1) those that measure time 
difference of sound travel between two points,°°4 and (2) those that measure the 
Doppler shift.°>& 


TIME DIFFERENCE 


Two transducers, alternately radiator and receiver, permit finding the dif- 
ference in transmission time for the two directions (fig. 1). Two radiators and two 
receivers allow pulses to be sent simultaneously in both directions (fig. 2); two 
closely spaced, parallel acoustic paths are used. 

A desirable feature of acoustic meters is their lack of moving parts. Out- 
put linearity is better than 1 percent and with this type is independent of the speed 
of sound in the medium assuming identical acoustic paths. With meters using only 
two transducers, this assumes the same acoustical path for a time duration greater 
than the double transmission time. In two parallel paths, this assumes identical 
conditions over the enclosed space during the simultaneous transmission time. 
Under conditions of strong thermal gradients, instantaneous velocity errors up to 
10 centimeters per second may occur.4 
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Wigure 1. Time-difference acoustic meter, one path. 
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Figure 2. Time-difference acoustic meter, two paths. 


Such temperature gradient effects would be serious if the meter were used 
to measure velocities associated with internal waves. Apparent velocity changes 
could be caused by the temperature gradients associated with the internal waves. 
The use of a meter with two transducers and a single path reduces this error. The 
typical ranges are from a few tenths of a centimeter per second to several hundred 
and overall accuracies 1 to 2 percent of full scale. However, there is a requirement 
for stable, expensive electronics because the differences in transit time are gener- 
ally less than a microsecond. The instrument is inherently capable of sensing 
direction as well as magnitude of flow, but its use in oscillating flow can be 
severely limited by flow disturbance. Vortex shedding from the transducers pro- 
duces high-frequency velocity fluctuations. Some meters may be zeroed in the 
field in the presence of flow, a necessary requirement for long-term measurements. 


DOPPLER SHIFT 


The Doppler acoustic meter depends on volume reverberation and a result- 
ing frequency shift (fig. 3). Acoustical energy is transmitted in a narrow beam 
from one transducer. A receiving transducer having a similar beam pattern is posi- 
tioned so that the two beams overlap. Typical volumes intersected are a few cubic 
centimeters at a distance of 50 centimeters from the transducers. The frequency 
of the reflected signal from the intersected volume filled with natural small scat- 
terers 1s compared with the transmitted frequency. The average velocity of the 
scatterers is assumed to be that of the fluid. The frequency difference is propor- 
tional to the fluid speed in the direction of the beam. 
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Figure 3. Doppler acoustic meter. 


Transmitted frequencies must be very stable; with typical values from 3 to. 
10 megacycles per second (Mc/s), the resulting Doppler shift is only a few hundred 
cycles per second (c/s). The method is capable of determining direction as well 
as magnitude. In view of the meter measuring the velocity some distance from it- 
self, little hydrodynamic disturbance might be expected, but this is only so when 
the meter is used in unidirectional flow and oriented to point into the flow. In 
oscillating flow, large distortions can result from disturbance of the flow about 
the case. Multiple sets of transducers mounted in a symmetrical case could con- 
ceivably be used, but would complicate matters considerably. 

Typical accuracies of the Doppler meter are 1 to 2 percent of full scale, 
and ranges are from a few tenths centimeters per second to several hundred. The 
computed flow rates are a function of the speed of sound, and corrections are 
needed. The directional response is unknown. The meter will respond to high- 
frequency fluctuations. 


SUMMARY 


The time-difference and Doppler-shift acoustic flowmeters are more suita- 
ble for constant-temperature, unidirectional flow, laminar or turbulent, than for the 
required need. Temperature effects (zero drift and multipath propagation); inability 
to measure and set the zero in some cases; hydrodynamic interference with the 
flow; and uncertain direction response pose serious problems. 


FLOW SPEED (KNOTS ) 


Impellor 


Impellor meters, used widely by industry, are employed by oceanographers 
in a configuration consisting of an axial rotor mounted on bearings inside a length 
of tubing with circular cross section. On the rotor are equally spaced blades 
placed at some angle to the pipe axis. The pipe gives the meter a degree of 
directionality. The angular velocity of the rotor is a function of the current 
velocity. This angular velocity of the rotor is measured indirectly by photoelec- 
tric or magnetic pickups mounted on the pipe to sense the passage of the blades. 
The number of pulses per unit time obtained from such pickups are converted to 
rotor velocity. 

The dynamic response for commercially available impellor meters is not 
given, but with low-mass rotors and low-entrained water mass, the response is 
sufficient for the required oceanographic conditions. ¢ 

A typical steady-state calibration curve furnished by Marine Advisers 
(fig. 4) is shown for their Model B-7C ducted current meter. Its specifications 
will be used for some calculations in order to determine if this type of meter is 
acceptable for the intended use. In order to depict the nature of the nonlinear 
output, figure 4 was replotted in figure 5, using linear abscissa and ordinate 
scales. In the region of interest (below 1 knot), the nonlinearity is due mainly 
to a threshold. Between the threshold velocity and 1 knot the output is nearly 
linear. 
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Figure 4. Calibration curve, Marine Advisers Inc. ducted meter model B-7C (five-bladed 
impeller); calibration date: 27 February 1963. 
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Figure 5. Marine Advisers Inc. ducted model B-7C (replotted from fig. 4). 
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MEAN VELOCITY 


Consider a hypothetical example of the effect of this nonlinearity when 
extracting a dc signal in the presence of an ac flow. Assume the ac flow is of 
low enough frequency so that the steady-state output calibration curve (figs. 4 
and 5) will hold. For velocities under 1 knot, let the output, e, be given by 


e=K(U—-U,) for U> Up > 0 (1) 
e=K(U+U,) for U < -Up (2) 
e=0 for -Up<U<Up (3) 
where 
kK = constant 


U = velocity 
Up = threshold velocity 
Let the actual velocity component, U, in the direction of the meter axis be 
U=U)+U,, sin wt (4) 


where U, = de velocity component in direction of meter axis 


U 


m 


amplitude of ac velocity component 


Now for the case where Uy > Up, the meter input and output appear as in 
figure 6. The mean, or filtered, output @, over one wave cycle can be found from 


A (GC PAAR 
call : 1 : 1 : 
e =o i) K(U - Up) d(wt) + = if K(U + Up) d(ct) + =| kU Up)d(at) (5) 

0 B D 
where 

Up - Up 
A=sin} = 
UF 
er Un + Ug 
B = sin'~ = 
Ore 
C=37-B 
D=3m1-A 
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Figure 6. Nonlinear distortion, impellor meter. 
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6=—— 
2 


[Uj(A+C-B+2n-D) +Up(C-B-A-2r+D) +U,,(cos B-cosC-cosA+cosD)] (6) 


To find the mean velocity, Ww, corresponding to the value of @, a formula 
similar to equation (1) would be used. Solving for w 


e+KU 
pails (7) 
K 


This would be the mean output of the real meter. It can then be compared with the 
velocity Uy. Although detailed curves are hard to plot because of the discontinuous 
nature of the output curve, a few values were calculated (Up = 4 cm/sec; Up from 
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5 to 15 cm/sec; and U,, from 5 to 20 cm/sec). Resulting errors of 20 percent in w 
values were typical. This effect is due to the filtering of the nonlinear output with 
a nonzero average and is not a function of frequency. Since only the static cali- 
bration is used, it is not an inertial effect. 

Steady-state linearity, threshold, and frequency response can be improved 
by various servo-controlled systems,®8 however, modifications of this kind would 
be difficult under laboratory conditions and probably impossible in the ocean 
environment. 


ROTOR ROTATION 


The impellor flowmeter usually cannot sense direction of rotor rotation and 
is therefore unable to sense current direction. The outputs from the pickups are 
usually the same regardless of the direction of rotor rotation. Several methods for 
the solution of the problem exist; one is given by Koppe.% 


ANGULAR RESPONSE 


Figure 7 shows the angular response for the Marine Advisers Inc. Model Q-8 
ducted meter, which will be considered typical for this kind of meter. The response 
is not a cosine function, and a single meter cannot be used to indicate the instan- 
taneous component of the real flow parallel to the principal axis without some error.* 

If more than one meter is used, the correct velocity components can still be 
found for some cases. Specifically, one possible case is where two-dimensional 
flow is present in the plane formed by the axes of two orthogonal meters. The de- 
tails, however, will not be given here since the more important case is where 
three-dimensional flow is present. Using three or more meters with noncosine 
responses does not appear to be feasible. 

The present need is to determine the seriousness of a noncosine response 
and to investigate the improvement of response characteristics. 


*Lane!0 states, for a meter of similar design, ‘‘. . . the proper component of flow along the 
axis of the instrument is recorded.’’ (Result of tow tank tests at constant speed.) However, 


no data were given and figure 7 will be assumed to be more typical for this kind of meter. 
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Figure 7. Directional response of Marine Advisers Inc. ducted current meter Q-8 
(data taken from Marine Advisers polar plot, operating instructions manual). 


SUMMARY 


The use of impellor meters for shallow-water oceanographic measurements 
is not recommended. Nonlinear outputs give rise to false de levels when filtered. 
Threshold and calibration changes at low velocities, caused by fouling, would 
probably be serious, and the noncosine response is a disadvantage. The use of 


two or three meters in the impellor system would result in hydrodynamic interference 


between the meters. 


Electromagnetic 


Electromagnetic flowmeters, or induction meters, have been towed from 
ships for oceanographic measurements or used to measure flow in channels. !! 
One possible oceanographic application concerns the flow in tubes with artifi- 
cially generated magnetic fields, where the tube gives the meter some directional- 
ity. A typical transverse magnetic-field configuration is shown in figure 8. Elec- 
trodes are mounted flush with the inside of the pipe wal! and are in contact with 
the fluid. The line joining the electrodes is perpendicular both to the magnetic 
field and the flow. Magnetic fields may be ac or de, although with electrolytic 
solutions such as sea water, ac fields reduce electrode polarization. 

Analytical work with these meters assumes that the fluid is nonmagnetic 
with the same permeability as a vacuum; conductivity of fluid, o, is assumed 
isotropic. The Hall effect is neglected, and variations in conductivity and 
thermoelectric effects are assumed negligible. The conductivity of the pipe walls 
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is assumed to be low compared with that of the fluid. Variations of electrical 
conductivity can generally be neglected if the impedance of the measuring instru- 
ment is high compared with the internal impedance of the meter. 


KLECTRODES 
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Figure 8. Transverse magnetic field configuration, electromagnetic flowmeter. 


DESIGN PARAMETERS 


The most important design parameters include choice of ac magnetic field 
frequency; length and diameter of pipe; choice of materials; effects of meter length 
onmagnetic field uniformity; methods of reducing and compensating unwanted ac 
pickup, and use of electrodes or coils for measuring induced voltage. One meter 
used in oceanography has a faired housing around the tube. 


ADVANTAGES 


The induction flowmeter used no moving or intricate parts. The method is 
inherently capable of measuring direction as well as magnitude of flow. Linearity 
is good; errors of less than one percent appear possible with careful design and 
proper use. The output is independent of the conductivity of the fluid if the mea- 
suring instrument impedance is high compared with the internal impedance of the 
meter. The output is therefore practically independent of changes in temperature, 
pressure, and density of the liquid. Without moving parts, upper frequency limits 
are high and the device can measure turbulent flow. In circular transverse field 
meters, concentric deposits of solid matter on the pipe walls, if they have the 
same conductivity as the fluid, do not change the calibration.!! The simple 
mechanical design allows for easy cleaning. 


DISADVANTAGES 


The induced voltage is usually independent of the electrical properties of 
the fluid. However, this is not true for all designs or conductivities.!2 If the 
meter requires salt water (conductivity, c= 4 mho/m) for accurate calibration 
instead of fresh water (tap water o = 10-2 mho/m), the choice of tow tanks is 
severely limited. Two or three orthogonally mounted meters might interfere 
hydrodynamically with the flow, depending on the size and shape of the sensors 
and on the distance between the electronics package and the sensors. 


DIRECTIONALITY 


A serious problem results from the use of the tube to obtain directionality. 
The direction response is unlikely to be a cosine function for two reasons: 


1. As the meter is turned through an angle © with reference to the flow, 
the mean volumetric flow will probably not be proportional to cos © for hydrody- 
namic reasons. 


2. As the meter is turned through the angle ©, the velocity distribution no 
longer remains symmetrical about the meter axis. 


Thurlemann (results summarized by Shercliff!+) proved that in a transverse - 
field flowmeter with circular cross section, the voltage induced in the electrodes 
is proportional to the mean volumetric rate of flow so long as the velocity distri- 
bution is symmetrical about the axis of the tube. Serious errors can result from 
flows with asymmetrical velocity distributions.8 Failure to produce a cosine 
response would lead to the same problems in resolving the instantaneous magni- 
tude and angle of the three-dimensional vector as was found with the ducted 
impellor meters. 


NONCOSINE DIRECTION RESPONSE. A noncosine direction response 
can represent a serious problem. Since instantaneous velocity is not required for 
this application, the question is posed as to whether the filtered output will 
indicate the steady flow correctly when given both an oscillating and steady flow. 
The noncosine response will be assumed to be the only fault of the meter. The 
meter will be considered perfect in other respects: linear output, perfect frequency 
response; and no zero drift. The ability to sense flow direction will be assumed. 


Let the meter be placed at a given depth with its axis parallel to the y-axis (fig. 9). 


A typical flow direction would be as follows: 


Let U be the instantaneous resultant flow at angle B to the y axis. The 
velocity components U,, U,, U,, and U, are shown. Let plane surface waves be 
traveling in a direction parallel to the U, component and give rise to the oscillating 
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Figure 9. Noncosine directional response, linear output example. 


part of the flow at the meter depth. Circular orbital velocities will be assumed. 
The flow components, including the steady-state values, will be given by 


U_= U, + U, cos wt (8) 
U,= Ug+ U,, sin wt (9) 
where U, and U, are the flow components shown, 
U, = steady-state flow component parallel to U, 
Uy = steady-state flow component in the positive z direction 
U,, = amplitude of the oscillating flow 
The output of the meter at angle 6 can be written as 
(10) 


e=e| a=o(R) 


where 


e| go = meter output at 6-0 


R = meter response at f 


OS Rall 
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Now since the output was assumed linear, equation (10) can be written 


O=IK (Ol Ine 

where 
Kk = constant (11) 
K>0. 


The mean (or filtered) value over one wave cycle, é@, can be found from 


21 
es IK 2 Vs : 
@=—— U R d(wt) where | U| = (U,“ + U*)? (12) 
on Z T 
0 


The value of R is a known empirical function of 8. The value of B is given by 


U, 
B = cos! ( 7 ) where U, = U, sin © (13) 


Using this value for U), equation (8) for U,, and equation (12) for U, 


(U, + U,, cos ct) sin © 


B=es! | (14) 
(U2 + U2)2 


As a check, let R = cos B. Using equations (12) and (14), where U,, U,,, and © 
are constants, 


1 51 | (Ui + Un cos wt) sin © 
a= || (2 Wa) Mra cae dwt = U, sin® (15) 
- (UzEU2)2 


Thus the correct steady-state component along the meter axis is indicated when 
the response is a cosine function. 

Equation (12) cannot be used at this time for an example of an electro- 
magnetic meter since an actual direction response is not known. However, to 
illustrate the method, a direction response the same as the Marine Advisers 
Model Q-8 ducted meter (fig. 7) will be used. This was a noncosine response 
that would probably be somewhat typical for an electromagnetic meter. The 
function of R vs § (fig. 7) was approximated by a series of straight lines, connect- 
ed by points taken at 5-degree intervals. Linear interpolation was used between 
points. The discontinuity at B = 76 degrees was avoided by letting R = 0.356 at 
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76 degrees and 0 at 77 degrees. In order to evaluate €, using this actual table of 
R vs B, equation (12) was modified for computer use as follows: 


N 


Eee 
e (16) 
NEA yy Ok 
n=0 


Substituting the value of U from equation (12), 


N 
e-— S. + (U2 + U2) (17) 
n=0 
where 
U, = Ug + U,, sin (nA) (18) 
U. = U; + Up, cos (nA) (19) 


R = function of 8 from the table 


{% _(Usin® 
B =cos (=) = cos (=—" ) 


© = positive-constant 


In equation (17), the positive sign is used if U) > 0 and the negative sign if U, < 0. 

The value of @ was computed for some typical values using A = 0.1 degree 
and N = 3600. This is the actual mean output of the meter. The ideal mean output 
along the meter axis would simply be U, sin ©. The results are plotted in figure 10. 
The ratio of the actual mean output to the ideal mean output is used as the ordinate 
and © as the abscissa. The asymmetric nature of the input and the noncosine meter 
meter response cause the errors in the mean output in the figure to range from 9 to 
23 percent. However, it should be noted from the curves that these errors are main- 
ly a function of © and almost independent of U; and Ug. This suggests the possi- 
bility of correcting errors of this type. A more detailed analysis should be defer- 
red until the exact nature of the response curve for an electromagnetic meter is 
determined. 
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EARTH’S MAGNETIC FIELD 


No tubes are used to obtain directionality in the method of induction-flow 
measurement that uses the earth's magnetic field. The induced voltages between 
electrodes can be given by 


V=UHL x<108 (20) 


where U = fluid velocity (em per second); H = earth’s magnetic field (gauss); 
L = distance between electrodes (cm); V = induced voltage (volts); and where 
U, H, and the line joining the electrodes are mutually perpendicular. 

The low strength of the earth’s field causes the signals to be very small. 
Large separation of electrodes thus becomes necessary to obtain reasonable 
voltage levels. A typical signal for electrodes of 100-meter separation is 2.5 
millivolts per knot. Since the earth’s field is stationary, the de voltages induced 
cause some polarization problems. 

Modern electrometers with very high input resistance are available to 
minimize loading problems. Calibration and testing in a tow tank is impossible 
for electrode spacings of this size. Location has to be distant from large bodies 
of iron such as tower structures or ships. Effects from the sea floor must be 
considered. Generally the velocity measured will be that component which is 
perpendicular to the line joining the electrodes and to the magnetic field. In 
practice, the line joining the electrodes can be placed parallel to the horizontal 
component of the earth’s field. Then only the vertical component of the earth’s 
field will produce an output. The output will be proportional to only the horizontal- 
velocity component perpendicular to the line joining the electrodes. Other orien- 
tations may produce signals due to vertical as well as horizontal-velocity compo- 
nents. The method is inherently capable of sensing direction of flow. The output 
polarity at the electrodes is merely reversed. Zero drifts can be of the same order 
of magnitude as the signal and are inherently difficult to determine. One method 
for shunting the electrodes to obtain an electrical zero is given by Mangelsdorf.!° 


MEASUREMENT OF VELOCITIES 


The measurement of velocities over large distances (100 meters) has a 
possibility for increasing signal-to-noise ratios. Assume the surface swell con- 
sists of plane waves of wavelength A. Let an instrument measure horizontal-flow 
components in a direction parallel to the surface wave propagation. Measure this 
flow over a distance of A or integral multiples thereof. Thus a rejection filter 
results for the major source of noise, the swell of wavelength A. The filtering is 
accomplished in the transducer itself and the method becomes Approach 3 (page 6 ). 
Typical values of A would be over 100 meters. In practice, both the direction and 
wavelength will change with time. Instrument placement would be determined by 
the expected swell direction and wavelength. 


SUMMARY 


It does not appear possible to place electrodes using the earth’s field to 
produce such a rejection filter. Perhaps another configuration employing ac fields 
could be constructed to give average velocity readings over large distances. 

The main advantages of the induction method are the lack of moving parts; 
high-frequency response; linear output; and direction-sensing ability. The major 
disadvantage with the tube-type meter using its own ac magnetic field is a probable 
noncosine direction response. The principal disadvantages with the earth’s field 
method are the low signal levels and zeroing problems. Installation and proper 
maintenance would also be a problem. Direction response is a cosine function, 
and the rejection-filter configuration should be kept in mind. 


Drag Force 


STEADY-STATE FLOW 
This method consists of recording the hydrodynamic drag on a fixed object, 
which for convenience is usually attached to a stiff suspension system. The small 
relative motion between the object and the mounting can be measured and recorded. 
The shape and size of the object and the suspension are design variables. 
For the steady-state flow in the region of interest here, the drag can be 
given by 
D =1/2 (Cp ApUlU) (21) 
where 
D = drag force 
Cp = drag coefficient 
A = cross-sectional area of object perpendicular to flow 
p = fluid density 
U = velocity of object relative to fluid 
The steady-state drag coefficient Cp for a given shape is determined 
experimentally as a function of the nondimensional Reynolds number. The Rey- 


nolds number, R, is given by 
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where 

U = velocity of fluid relative to object 

d = characteristic dimension of object 

u = dynamic viscosity 

v = kinematic viscosity 

p = fluid density 

As the Reynolds number increases, vortexes are produced and shed behind 

blunt bodies such as cylinders. When R is between 10 and 40, stationary vortexes 
are formed. Stable and well-defined vortex patterns are carried away or shed for R 
between 40 and 150. The line of separation moves back and forth, and transverse 
forces appear which change direction with time. Using the Strouhal number, which 


is nondimensional and determined experimentally, the frequency can be found by 
the formula! 


= Sill) 
——— 9 
f 7 (23) 
where 
d = diameter of body 


f = frequency of oscillation 


U 
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velocity 
S = Strouhal number 


As R is increased from 150 to 104, the shedding becomes irregular. 


SHALLOW-WATER APPLICATION 


Considered so far have been objects in a steady-fluid flow, or objects 
moving at fixed velocities through a still fluid, after equilibrium has been reached. 
In the shallow-water case, however, oscillating flows are present because of wave 
action. This motion means that fluid acceleration is present which also influences 
the coefficient of drag,!5-!7 the wake formation,!® and the apparent mass increase, 
which is called virtual, or added mass.!9-22 

Of prime importance are the variations in drag coefficient and mass coeffi- 
cient with fluid acceleration. This problem has been studied theoretically and 
experimentally for unidirectional accelerated motions, arbitrary motions, and 
sinusoidal oscillations. Keulegan and Carpenter!® give a good review of these 
references, and present results of sinusoidal fluid oscillations for cylinders and 


plates. Their work is a good reference because it shows a systematic correlation 
of the variation in the coefficients of drag (Cp) and mass (C,;). The following 
items are described: analysis of forces on objects at rest in a moving fluid, with 
the special case of a cylinder in a sinusoidal fluid velocity; the experimental 

wave tank and force dynamometer; data reduction; flow patterns which resulted; and 
possible reasons for the variation in coefficients due to vortex formation. The 
results of primary interest were for cylinders: 


1. The average values of Cp and Cy, during a wave cycle were correlated 


m 


with the period parameter where U,, is the amplitude of the sinusoidal fluid 


velocity. TT = period, d = cylinder diameter. No correlation was found with the 
Reynolds number. The results are plotted in figures 11 and 12. 
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Figure 11. Variation of drag coefficient of cylinders. 15 
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Figure 12. Variation of inertia coefficient of cylinders. 15 


2. The variation of C\, and Cp throughout a given wave cycle were 
experimentally determined for use in the following equation: 


ae ae : 
ogp t 1/2 (Cp d pulU) (24) 


= Cy pA 
where 
F = force per unit length in direction of flow 


C\, = coefficient of mass 


po = fluid density 
2 


Ag = — where d = diameter of object 
U = undisturbed fluid velocity 
Cp = drag coefficient 


See results in figures 13, 14, and 15 for different values of period parameter. 
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Figure 13. Example of variation of the inertia and drag coefficients of a cylinder 
during a wave cycle. 
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Figure 14. Example of variation of the inertia and drag coefficients of a cylinder 
during a wave cycle. 
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Figure 15. Example of variation of the inertia and drag coefficients of a cylinder 
during a wave cyclel 


3. The maximum force and phase were found as a function of period 
parameters (fig. 16). 


The problem with using equation (24) for an actual current meter is that 


] 


both U and are unknown. It is therefore highly desirable for the parameters to 


be chosen so that the acceleration term can be neglected. Furthermore the value 
of Cp should be constant over the operating range of the instrument. Corrections 
for variations of Cp would be unrepeatable since these variations are connected 
with fluid acceleration and eddy formations, and varying coefficients complicate 
data reduction. Steady-state data show that the range of Reynolds number for 
constant Cp is between 1 x ae 2105. Figure 11 shows that the higher 
values of the period parameter ee ) are desired because Cp and C,, are nearly 
constant in this range. In the shallow ocean, U,, and T are fixed by the amplitude 
and period of the swell. The only variable over which there is control is the 
diameter d of the cylinder, which must be small if period parameters are to be 
large. The reduction of d, however, reduces the drag force proportionately. These 
smaller forces produce more demands on the electrical transducer. 
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Figure 16. Variations of the magnitude and phase of the maximum force 


on eylinders. 15 
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DESIGN PARAMETERS 


An actual meter of this type (fig. 17) was constructed in order to learn 
more about the range of variables to be encountered. Data samples (figs. 18 and 
19) were taken at the NEL Oceanographic Research Tower2? where the water 
depth is 60 feet. Both records were taken at a depth of 40 feet; one shows the 
north-south component and the other the east-west component. These records 
will serve as examples of typical parameters that influence the accuracy of the 
method. Keulegan’s results can then be applied. 

A cylindrically shaped object attached to a cantilevered beam (fig. 20) 
was used. The beam shape causes the cylinder to bend in only one plane and 
therefore the meter is directional. The bending is measured by strain gages 
cemented to each side of the beam. The hydrodynamic drag on the cylinder 
causes the beam to bend and compress one strain gage and stretch the other. 
This produces an equal resistance change in the gages, which in turn produces 
a voltage output in a de Wheatstone bridge circuit. 

The problems with the instrument include how to make a sensitive 
electrical transducer with satisfactory, long-term stability and a low temperature 
coefficient while retaining rugged mechanical construction and a high natural 
frequency. However, consider now only the inherent hydrodynamic characteristics 
and assume that a perfect electrical transducer can be built. Its voltage output 
would be proportional to the drag force on the cylinder and would have perfect 
stability without temperature coefficient. 

el 


First the values of period parameter 77 are computed. The approximate 


range of period parameter during the 3 days when data were taken was 40 to 240 
for the north-south direction and 60 to 410 for the east-west direction, with d 
constant and equal to 0.51 centimeter. The ranges cover about 90 percent of the 
cases. A typical value for the period parameter was about 200 with T = 10 seconds 
and U,, = 10 centimeters per second. 

The easiest and most accurate quantity to find is the maximum velocity 


] 
during a wave cycle since at this point = = 0. Figure 16, which shows the 


variation of the magnitude and phase of the maximum force as a function of period 
parameter, can be used. These curves assume a sinusoidal velocity variation. 
For example by using a typical value of 200 for the period parameter, 

a : 


_'m__ approaches 0.6. Therefore F = 0.6 pU,,2d, which is the same as using 


fe) ime d 
only the second term in equation (24) with a drag coefficient value of 1.2. This 
is the steady-state drag coefficient value commonly used. Thus for high period 
parameter values, the acceleration term is relatively unimportant in finding 
maximum velocity. Note that the phase shift of the maximum force does not 
approach zero even for large-period parameter values. 

Consider the case where it is desired to find the velocity throughout one 
wave cycle. As shown in figures 13 and 15, if the period parameter is either 
small (3) or large (44.7), the values of C,, and Cp are roughly constant over the 


het ; t 2 = 
entire cycle. The greatest variation of C,, and Cp is at 7 = 0.25 and 0.75. 
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Figure 17. Drag force current meter. 
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Figure 18. Recorder sample during run 2 (north-south). 
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Figure 20: Schematic drawing of beam and cylinder. 
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These average values (for one cycle) can be found from figures 11 and 12 for 
use in equation (24). 
: du f : seas 
In practice, both U antl are unknown, and knowing Cy, is of little use. 


The errors in neglecting the acceleration term can be large, even at large values 
Ua 
of period parameter. Take the typical case of ee 200, T = 10 seconds, 


U,, = 10 centimeters per second, d = 0.51 centimeter. Extrapolating figures 11 
and 12, Cp) = 1.2, and Cy, = 2.6. Substituting these values into equation (24), 


dU 
[n= aa + 0.31 U|U| (25) 


Let U=sin wt. Values of F from the above equation are plotted in figure 21, 
first using only the velocity term and then using both terms. The force curve 
using both terms is what a real current meter would actually detect. Therefore 
by using this force curve and solving for U in equation (24), where the accelera- 
tion term is neglected, the indicated velocity can be plotted (fig. 22). The true 
velocity, U = sin wt, is also plotted for comparison. The large errors are obvious 
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from the curves. Making larger, by making d smaller, probably will not 


change the values of C,, and Cp. This is a reasonable assumption for Cp until 
the critical Reynolds number is reached. The assumption for Cy, may not be valid 
although the curve in figure 12 appears asymptotic for large-period parameter 
values. The ratio of the velocity term to the acceleration term varies as 

1 ; 

qin equation (24). By making d smaller, the acceleration term could be made 
less significant. Making d several times smaller, less than 0.51 centimeter, would 
reduce the forces to such a level that they would be hard to measure. 
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Figure 21. Drag force current meter example. 
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Figure 22. Drag force current meter example. 


For values of period parameter near 15 (fig. 14), neither C,, nor Cp can 
really be considered constant throughout the wave cycle, and equation (24) 
cannot be used with constant coefficients. If data had period parameter values 
near 15, additional experimental work would be necessary. 

An effort to recover any low-amplitude, low-frequency, or essentially de 
signals in the presence of these typical swell records places very severe require- 
ments on a system. No de velocity signals were ever detected visually on any 
of the Tower records. 


ACCURACY 


Consider now the accuracy of finding the steady d-c velocity after 
filtering out the ac terms. For the first case, assume that high-period parameter 
values are used and therefore the drag coefficient will be constant throughout 
the wave cycle. If the inertial term in equation (24) could be neglected, the 
actual force per unit length on the cylinder, fF’, would be given by the second 
term of equation (24), 


F = Ky U|U| (26) 


where Ky = 1/2 (Cp d p) = constant. 


Let the actual velocity U = Up + U,, sin © where Up and U,, are constants (fig. 23). 
Since the meter output, e, can be made linear with F, 


e = Ky Kg U|U| (27) 
where K3 = constant. 
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Figure 23. Ac + de flow throughout one wave cycle where inertial term may be 
neglected and drag coefficient is constant; drag force current meter. 


Now this actual output could be converted to the corresponding U for every value 
of © before any filtering is done. Thus the indicated velocity, U;, would simply be 


72 . 
= +1fU>0 : 
eae e (28) 
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Up=— | Upde=Up (29) 
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Therefore, for this case, the correct mean velocity Ug would be found. 
Next consider the case where the inertia term is still assumed negligible 
when solving for U; (when in fact it is not). The drag coefficient is still 
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considered constant. For this case the actual force must be given by both 
terms of equation (24), 


5 ll) 


where 
K, = 1/4 (Cy om d?) = constant 


As before, let the true velocity be given by U = Up + U,, sin 9 and the actual 
output be linear with force. 


: dU 
Thus e=KykK3 as KoK3 U|U| (31) 


Since the inertia term is still assumed negligible, equation (28) is again used to 
solve for the indicated velocity, resulting in 


Wy 
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And U; can be found, 
21 
= eg x 
Nt = i) U;,d9 (33) 
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This method will not normally give the correct mean value Ug when Up 4 0. As 
examples, values of U, as found by numerical integration for two cases of Ug/U,, 


are shown below. 

The values of kK, = 0.54 and Ky = 0.31, obtained earlier for the typical 
example U,, = 10 centimeters per second, T = 10 seconds, d = 0.51 centimeter, 
were used. 


Percent Error 
Between 
Up and U, 


The indicated values of U, were lower than the correct value Ug, and quite sizable 
errors resulted. The values through the wave cycle for the case where Up/U,, = 0.2 
are plotted in figure 24. 

Other shapes can be used instead of cylinders; however, the same type 
of problem still exists. The inertial term cannot be neglected for flat plates except 
at very small values of period parameter.!® Adding to this is that quite large changes 
in Cp and Cy, throughout the wave cycle exist for all values of period parameter. 
Although spheres were not studied,!5 the results would probably be similar to those 
with cylinders. 
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Figure 24. Drag force meter example. 
ADVANTAGES 


The good features of the drag force meter from a hydrodynamic standpoint 
are as follows: 


1. The meter is small and has few moving parts. 
2. The directional response can be made a cosine function. 


3. The output will indicate both the magnitude and direction of a given 
velocity component. 
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DISADVANTAGES 


The bad features of the drag force meter froma hydrodynamic standpoint 
are as follows: 


1. The meter has a nonlinear output, even in steady flow. 

2. Cp and Cy, vary for oscillating flows. 

3. The values of Cy, are not small compared with the values of Cp with 
dU/dt unknown for data reduction use of inertia term. The neglect of this term 


may result in sizable errors in solving for the dc level. 


4. The slope of the force (or output) vs velocity curve is zero at zero 
velocity, which makes the zero inherently difficult to determine accurately. 


SUMMARY 


The drag force method is not suitable for measuring low-level dc signals 
in the presence of these typical swell conditions. Even if a perfect electrical 
transducer could be built, the inherent hydrodynamic problems would be too great. 
This method is, however, being used successfully as a ‘‘surge recorder’ where 
only the maximum value of the oscillating flow is required. 


LAGRANGIAN METHOD 


The Lagrangian Method involves determining the position of artificial 
tracers in space and time. Von Arx? lists, as types of tracers, drift bottles, 
drift cards, radio buoys, contaminants such as dye or radioactive waste, current 
poles, deep drogues, neutrally buoyant floats, and ship drift. Current poles and 
‘*deep’’ drogues can be used in shallow water. Even though measurements are 
taken in shallow water, it is still desired to measure velocities at some predeter- 
mined depth. Near the surface, current poles could possibly be used whereas deep 
drogues would be needed for other depths. Current poles are basically long, thin, 
vertical spar buoys, with only a small portion above the water line, which are 
used for visual tracking. In the deep-drogue configuration, a surface float is 
connected to a submerged object having large drag by a fine line. The float is 
tracked visually. 

Of the serious faults and disadvantages inherent in the Lagrangian Method 
when used in shallow water, the most obvious is the necessary support effort. 
Setting, tracking, and picking up drogues is by no means an easy task either at 
night or for long periods of time. 

The output, which is a function of space and time, is not electrical and 
might not comply with the condition of a point measurement satisfactorily. Since 


the output is sampled, low-frequency aliasing must be kept in mind. Optical 
tracking consists of sampling positions, usually at equal time intervals. 

The drogue, which must respond to the horizontal components of all 
velocities present, should either follow the de and ac velocities perfectly or act 
like a low-pass filter— have a reduced response to the ac flow, but follow the de 
flow perfectly. However, the drogue would not normally follow the ac flow 
perfectly. A difference in velocity between the drogue and the ac flow will exist 
and vary with time. The drogue is a passive device that requires an external force 
obtained from this velocity difference in order to operate. The resulting force, 
however, on the drogue is not linear with velocity and can generate a false de 
level except when the true de level is zero. With the fixed drag-force meters 
equipped with electrical readout, the nonlinear force vs velocity effect could in 
principle be taken out in the process of data reduction, but this is impossible with 
drogues. 

Assume a velocity difference in the horizontal plane, U, between the drogue 
and the water given by U = Ug + U,, sin ©, where Up and U,, are constant. Assume 
that equation 24 can be used to describe the resulting force on the drogue. Assume 
Kk, and Kg are constant throughout the wave cycle (this might not be a valid as- 
sumption in actual practice). Also assume that U,, > Ug. The net force per cycle, 
F, can be found. 
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The mean value, F, was evaluated numerically for one set of values. For a ratio 
of U,,/Up = 5, the net force per cycle, F, was about 6.4 times larger than the ideal 
force which was found by using only Up. False de levels can thus be generated 
by filtered nonlinear ac terms and was in this case a function of Un UG abhis 
was only a crude example to show that the effect can occur. A more complete 
description would be quite complicated. Model studies would probably be 
necessary to show how serious the effect would be with real drogue configurations. 
Another disadvantage appears with drogues constrained to stay at one 
preset depth. At any given time and at a constant depth, some horizontal velocity 
distribution exists. Assume that this is a sine wave. This sinusoidal spatial 
distribution travels through space as a function of time at the phase velocity, c, 
of the surface wave. Since the drogue is free to move only in the horizontal plane, 
the condition for the net horizontal force on the drogue to be minimum would exist 
only when the drogue is traveling at the velocity, c, also. This condition would be 
independent of the initial depth of the drogue as long as it stays at this depth. 
The use of drogues in shallow water is not recommended. 
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Short-Distance Tracers 


Normally the Lagrangian Method is used to describe the tracking of artifi- 
cial tracers over relatively large distances (e.g. miles). However, several low- 
inertia tracer methods exist for shorter distances (e.g. less than one meter). 
Three methods are briefly listed: 


1. Nuclear magnetic resonance .2° 


2. Photographing the motion of small drops of colored liquid. 2° 


3. Measure the temperature fluctuations in the wake of a sinusoidally 
heated fine wire.2/ 


These methods would be more suitable for laboratory or industrial appli- 
cation (e.g. pipe flow) than for oceanographic field work. Problems in the design 
of a three-dimensional ocean meter would be large and mainly practical in nature. 
These methods are not recommended at this time, but should be kept for possible 
future use. 


MISCELLANEOUS 


Savonius Rotor Meters 


These should never be used in oscillating flow because of poor dynamic 
response and rectifying action due to the omnidirectional response.29 


Thermo-Anemometer Meters 


These types, such as the hot-wire anemometer, relate fluid velocity and 
the quantity of heat lost by a heated measuring element (thermistors, wires, films). 
They are inherently suited for high-frequency velocity fluctuations and not for 
applications requiring long-term stability. Changes in sensitivity due to fouling 
would probably be serious. Thermo-anemometer meters have no inherent capability 
for sensing direction and would require shrouds and/or multiple sensor arrays. 
They are not recommended for the present application. 


CONCLUSIONS 


1. The acoustic flowmeters, for several reasons, are not suitable for 
measurement of steady-state currents in shallow water. 


2. The nonlinear outputs and the noncosine response of impellor 
flowmeters give rise to false dc levels when filtered. 


3. The inherent hydrodynamic problems of the drag force method are too 
great. 


4. The necessary support effort required for the Lagrangian method makes 
its employment impractical. Nonlinear forces produce false de levels. 


5. The electromagnetic flowmeter is inherently capable of measuring 
direction of current; the output is linear. Direction response is unknown for the 
tube configuration. 


RECOMMENDATIONS 


The electromagnetic method is recommended for the intended use. 
Specific recommendations are as follows: 


1. Investigate new configurations using ac fields that will provide 
minimum hydrodynamic interference and possibly measure flow over large dis- 
tances. 


2. AC Magnetic Field Tube Configuration (fig. 8): Experimentally deter- 
mine the directional response of this type. Investigate effects of different 
length-to-diameter pipe ratios and pipe fairing shapes in order to obtain a nearly 
cosine response. By using the best response curve obtained, calculate the 
expected errors (similar to example on p. 25). Determine if these errors are 
acceptable. 


3. Earth’s Field Type: Use of this method is recommended only if a 


feature for zeroing the meter is employed (similar to description by Mangelsdorf!2). 
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